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Corrections to the universal results 
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We investigate the deviation of the level-correlation functions from the universal form for the 
chiral symmetric classes. Using the supersymmetric nonlinear sigma model we formulate the per- 
turbation theory. The large energy behavior is compared with the result of the diagrammatic 
perturbation theory. We have the diffuson and cooperon contributions even in the average density 
of states. For the unitary and orthogonal classes we get the small energy behavior that suggests a 
weakening of the level repulsion. For the symplectic case we get a result with opposite tendency. 
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I. INTRODUCTION 

Statistical properties of disordered systems have at- 
tracted interest since the work of Wigner and Dyson. 
At present, we recognize that disordered systems_are 
classified with the notion of the symmetric space.u In 
addition to the three traditional classes of disordered 
systems — unitary, orthogonal and symplectic — we have 
three classes for chiral symmetric systems and four classes 
for normal-superconducting systems.H 

The chiral symmetric Hamiltonian can be written as 

tt _ ( uji - iiv 2 \ m 
H ~[u 1 +iu; 2 J' (1) 

where uj\ and a->2 are random Hermitian matrices. This 
Hamiltonian is relevant for the motion of a -single elec- 
tron in a lattice with random magnetic fieldsjj which has 
been under intensive study because of its relation to the 
fractional quantum Hall effect close to half filliiagj and to 
the gauge theory of high T c superconductivity.^! Further, 
the low energy properties of QCD are also believed to be 
described by this type of Hamiltonian.El 

This Hamiltonian has pairs of eigenvalues with ±E. 
The zero energy point-becomes a special point and new 
universality appearsBQ This universal nature at the er- 
godic limit is studied very well and we know that even 
the average density of states becomes universale 

In the ergodic limit, we are interested in phenomena 
whose time scale is much larger than the diffusion time 
1 /E c where E c is the Thouless energy. Hence an electron 
has enough time to wander everywhere in the system and 
the spatial structure of the system can be ignored. If 
we want to discuss relatively short-time phenomena (or 
long-range correlations of energies), we need to leave the 
ergodic region and enter the diffusive region. 



In the diffusive region of the traditional three classes 
of disordered systems, it is well known that the diffuson 
and cooperon modes that express the quantum interfer- 
ence effects play an important role. Especially owing to 
these modes there appear the deviations from the uni- 
versal behavior. The aim of this paper is to give analo- 
gous results for the chiral symmetric disordered systems. 
The deviations from the universal form are interesting 
from the viewpoint of the weak localization because we 
do not know whether or not the localization occurs for 
these classes. Besides, mesoscopk. devices with random 
magnetic fields are now available .□ 

In our previous work,E2l we discussed these problems 
in the context of QCD for the chiral unitary class but 
some of the result contain some errors. In this paper, for 
all three - chiral orthogonal, unitary and symplectic - 
classes, we calculate the level correlation-Junctions. Us- 
ing the Efetov's supersymmetry methodjii we formulate 
the perturbative expansion and calculate the large energy 
asymptotics of the level-correlation functions. We derive 
the small energy behavior of the density of states us- 
ing theJmproved calculation developed by Kravtsov and 
Mirlin.EZl The large energy behavior is compared with 
the results of the diagrammatic perturbation theory to 
examine which contributions are important and which 
ones make differences between the chiral classes and the 
traditional ones. 



II. LEVEL CORRELATION FUNCTIONS 

We formulate the calculation of the level-correlation 
functions by using the supersymmetry method. We con- 
sider the effective action 

S = ^KvJ ^ tr gQW^ + ^7 / ^trglVQ^)] 2 , 
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where trg denotes the supertrace0 and D is the diffu- 
sion constant and A is the mean level spacing. V = L d 
is the volume of the system and T, z = diag(l, 1,-1,-1). 
We note that this is the effective action for the one-point 
function. The supermatrix Q(x) that parametrizes the 
saddle-point manifold is a 4 x 4 supermatrix £or the chi- 
ral unitary class and satisfies Q 2 (x) = — 1£3 For the 
chiral orthogonal and symplectic classes, the size of the 
matrix is duplicated. The structure of the superma- 
trices comes from chiral symmetry and supersymmctry 
(and the complex conjugation property). This nonlinear 
sigma model is-derived from the schematic model of ran- 
dom matrices. E£ -A. similar model is derived by using the 
replica method.li3 Related models using the supersym^ 
metry method are derived for the random-flux model,EI 
for a Dirac particle in gauge field disorderE.3 and for a 
two-dipensional(2D) electron gas in a random magnetic 
fieldJij Although Eq.(||) is derived starting from a partic- 
ular model, we think this form is broadly valid because we 
know, from the experience for the traditional classes, the 
effective action does not depend on the details of specific 
models. In the case of d = 4, a similar model has-, been 
discussed in great detail in the context of QCDE3 The 
nonlinear sigma model was written down based on the 
symmetries of the theory and the nature of the diffusive 
regime as well as the ergodic regime was studied. 

The average density of states can be calculated from 
the following expression: 



(P(E)) = 



4AV 



Im / DQ 



d d xtrgkQ(x)Y, z 



(3) 



where k — diag(l, —1) ® I2. For the ergodic limit (which 
means we ignore the spatial dependence in the above ex- 
pression), by parametrizing the saddle-point manifold, 
we can calculate the density of states and get the uni- 
versal result.E3 In this paper, for the formulation of the 
perturbative expansion, we take another parametrization 



Q(x) = -iE,(l + P)(l 

= -i S 2 [ 1 + 2P + 2P 

it{x) 
-it(x) 



py 1 

2 



2P 6 + 2P 4 



P 



(4) 
(5) 



where t(x) is a 2 x 2 (4x4) supermatrix for the chi- 
ral unitary (orthogonal and symplectic) class. This 
parametrization has the advantage that the measure is 
normalized. A similar parametrization— is used for the 
traditional classes xl In previous works ,113 we have taken 
another parametrization. The measure is not normalized 
and has been treated perturbatively. 

Using the parametrizaton (^), we can formulate the 
perturbation theory. The free action can be written as 



SQ = SrJ ddxtIgp2 ~ iv/ ddxt ^pf. (6) 

The perturbative expansion is expressed using the diffu- 
sion propagator 



IL(E;x,y) =J2^(E,k)e ik ^, 

k 

1 



U(E,k) = 



(7) 
(8) 



A. Large- E behavior of the density of states 

The level-correlation function can be expressed in pow- 
ers of the diffusion propagator. This expansion is valid 
for the large E (large E/A). Therefore, we can see the 
large- E asymptotics of the level-correlation function in 
this calculation. The calculation goes along the same 
lines as the traditional class. The only difference is the 
contraction rules. For the chiral unitary class, the con- 
traction rules are expressed as 



{trg AP(x)BP(y)) l 



1 



IL(E;x,y) 



trg A trg B 



(trg AP(x) trg BP(y)) ( 



-trgAS z trgSS, 
-trg trg BE, 

+trgAE y trg£E y 



= ±IL(E;x,y)trg[AB 



(9) 



(10) 



where A and B are arbitrary 4x4 supermatrices and 
(• • -)q denotes the average with respect to the free ac- 
tion Sq. S x , E y , and T, z are 4x4 Pauli matrices. 

We can see-.that the differences from the traditional 
unitary classtll are the third and fourth terms of Eqs. 
@ and (|l0j). These differences come from the structure 
of the matrix P. For the traditional unitary class, the 
matrix P can be written as 





-it 2 i(x) 



it 12 (x) 




(11) 



kt 12 . 



The matrices ti2 and t 2 i relate each other with tzi 
When we take the contraction, these matrices couple 
with each other but there are no self-couplings such as 
(Ati 2 Bti2) . Due to this difference, we get the above re- 
sults for the contraction rules. 

Using these rules, we obtain the density of states for 
the chiral unitary class as 



(p(E)) =±Re 



1 



2{£n(£,fc)} 2 



2 



16wiE 
A 



(12) 



For the ordinary class, this function just equals to 1/A. 
The nontrivial contribution appears due to the chiral 
structure. We note that the first-ordex_correction (II 1 
term) vanishes. In the previous work ,113 we have erro- 
neously neglected the higher-order terms coming from 
(VQ) 2 and have got a wrong result. A similar term in 
the present parametrization gives no contribution. 
If we take the k = mode only, we have 



(p(E)) = £ 



1 



8(TrE/Ay 



(13) 



which coincides with, an asymptotics of the exact result 
at the ergodic limit E 

Similar calculations can be performed for the chiral 
orthogonal and symplectic classes. For the chiral orthog- 
onal class, the saddle-point manifold is parametrized by 
the 8x8 supermatrix Q(x) . We have additional condition 



Q = CQ C 7 



(14) 



where C is the complex-conjugation matrix. Similar cal- 
culations apply for the chiral symplectic class byrthe re- 
definition of the complex-conjugation matrix C.til The 
contraction rules can be written as 



(tig AP(x)BP(y)) Q 



-U{E;x,y) [trg A trg P 

-trgAS^trgSS, 
-tr gJ 4E x trg.BE, 
+txgAE y trgPE y 
+a trg (AT, X CB T C T T, X 
+AY, y CB T C T Y, y 



+ay, z cb t c t y: 



-ACB T C T * 



(15) 



i trg AP(x) trg BP{y)) Q = ~ U(E; x, y) trg 



AB 



+AY> y BY>y - ACB T C T 

+AZ Z CB T C T Y, Z 

+AY, X CB T C T Y, X 

(16) 



where a = +1(— 1) for the chiral orthogonal (symplec- 
tic) class. The presence of the terms including E, and 
E y is due to chiral symmetry. The terms including the 
complex conjugation matrix can be interpreted as the 
cooperon contributions. Using these rules, we obtain 



(p{E)) = ~Be 



l + 2a ^2fl(E,k) 



(17) 



The first-order contribution appears in contrast with the 
unitary case. 

If we take the zero mode only, we have 



(p(E)) =iRe 



2irE/A 



(18) 



We see that the 1/E term does not contribute to the 
density of states because this term is imaginary. On the 
other hand, this is not the case if one includes the k ^ 
modes. The second term in Eq.(|l7j) survives and repre- 
sents the lowest-order quantum-interference effect. We 
can see the signs in the two classes are opposite, which 
is reminiscent of the weak localization effects of the con- 
ductance for the traditional classes. The diffusion mode 
is essential for this effect to appear. 



B. Small-iJ behavior of the density of states 

The above results can not describe the small- E behav- 
ior because we do not treat the zero mode properly. For 
the small energy, the diffusion propagator diverges at the 
zero momentum, which means the breakdown of the per- 
turbation. This defect can be-d|Kjumvented by consider- 
ing the zero mode separatelyllJ'ta By this improvement, 
we can discuss the small- E behavior of the level correla- 
tion functions that is highly nonperturbative. 

The parametrization of Q is as follows; 



Q(x) = T^Q(x)T . 



(19) 



The zero mode Tq is treated exactlyS and the nonzero 
modes Q(x) are treated perturbatively as before. 

After the perturbative calculation of the nonzero 
modes, the average density of states is expressed as 



(P(E)) 



1 

A 



Re / dQo exp 



ixiE 



trg^o 2 



l+-trgfcP 2 



A 

i + 2n 2 (o) 



2niE „, . . o 
II 2 (0)trgP 2 

AniE ( 1 
V 



A 

2w 2 E 2 ( 1 



A 2 V V 



J d d xU 2 (x)^j trgP 2 
J d d xU 2 (x)^j (trgP 2 



(20) 



where H(x) is the diffusion propagator with E = and 
Po is the supermatrix which parametrizes the zero mode. 
This expansion is valid for E <C E c = D/L 2 because 
Eq.(p0[) is expanded in powers of E. 

The zero mode integration is performed and we obtain 
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- (P(E)) 

7T 



,(1)2 



' (1)2 (2)' 
a d , a d 



8g 2 Ag 2 / dz 



0(l/g 3 ) 



Ps(z) 
(1)2 



J 2 (z) + J 2 (z) + ^ZJ 2 (Z) 



(2) 



where z = nE/A, g — -kE c / A and p s (z) is the univer- 
sal function at the ergodic limit.i Constant a^ 1 ' 2 "* is the 
momentum integrations of the diffusion propagator 



-y 

n=£0 

-y 



1 



n#0 



(n\ + 



(22) 
(23) 



These constants depend on the spatial dimension d 
a[ 1} = 1/6, 4 X) = i InL/Z, - L/Z, af 5 = 1/90 z 



0.0111, a 



(2) 



0.0266, <4 2) 



0.0527 (/ is the mean free 



path that decides the cutoff scale). For some cases, e.g., 
the renormalizatoin group calculation, it is more con- 
venient to include counterterms of higher order in the 
action (||) instead of using cutoff parameters explicitly. 
For the small z (z <C 1 <C g), we have 



7T A 



,(1)2 



,(2) 



4 4 g2 



For 1 « z « .g, 



A(p(E)) 



,(1)2 



,(2) 



(25) 



Eq.(p5[) can also be obtained from Eq. (|l2|) . 

We can perform the same calculations for the other 
classes. In this case, there remains a 1/j term and we 
get 



- (P(E)) 



Ps{z). 

(26) 



It is remarkable that for all three classes the deviations 
can be expressed with the universal function p s (z) and 
its derivative. Whether this also holds for higher-order 
terms seems an interesting problem. The universal func- 
tions p s (z) for the orthogonal and symplectic cases ateal-. 
ready derived by the orthogonal polynomial methodO'EHI 
For the small z, we have 



A 1 / a (1) \ 

- (P(E)) ~ g ( 1 + -fg~ J ' ortho S onal > 



(P(E)) 



For 1 « z « g, 



2a 



(i) 



A(p(E)) 



z , symplectic. 



4 5 ' 



(27) 
(28) 

(29) 



(21) 

which is again consistent with Eq. 

Equations (24) and (|27j), and ( |28| ) mean that the re- 
sulting correction does not change the power behavior of 
the density of states, but renormalizes the correspond- 
ing prefactor. Similar expressions exist for the tworfioint 
level-correlation function for the traditional classes hi We 
obtain the results for three traditional classes that are 
interpreted as a weakening of the level repulsion. The re- 
sults for the chiral unitary and orthogonal classes could 
be interpreted similarly as a weakening of the repulsion 
between a pair of levels with the energy ±E. The re- 
sult for the symplectic case does not accept this trivial 
interpretation and is interesting. 

Gade analyzed a similar model in the thermodynamic 
limit by the renormalization group method for the chiml 
unitary class and got the divergence at the band centertj 
for d = 2. Our analysis implies Gade's result. It is inter- 
esting to perform the renormalization group analysis for 
the chiral symplectic class since our result suggests the 
vanishing at the band center. 

We note the dimensional dependence of the conduc- 
tance g. For A ~ L~ d and E c ~ L~ 2 , we have 
(24) .9 ~ L d ~ 2 . Therefore, our 1/g corrections become ir- 
relevant for d > 2 in the thermodynamic limit. 



C. 2pt. level correlation function 

We calculate the large-E behavior of the two-point 
level-correlation function. In this case, the saddle-point 
manifold is parametrized as 



P = 





-it(x) 



Q(x) = -i£*(l + P)(l - P)- 1 

it(x)\ s _ / h{x) t 12 (x) 
) ' t[X) ~ Ui(aO t 2 (x) 



(30) 
(31) 



Here, 2x2 supermatrices t\ and t 2 are due to the chiral 
structure of the model and t\ 2 and . t<u fire due to the 
calculation of the two-point function .113113 

For the unitary class, the leading-order contribution of 
the connected part of the level correlation function comes 
from the following contraction:E£l 



tTgkt 12 t 21 (x)trgkt 21 t 12 (y) ) = 4U 2 (E;x,y). (32) 

Q 



We have 



4 



A 2 [(p(E l )p(E 2 ))-(p(E l ))(p{E 2 )) 
Y / [ti 2 (E,k)+U 2 (-E,k) 

k 

+n 2 (w/2,fc) + n 2 (-w/2,fc) 



(33) 



We see that the last two terms already appear in the 
ordinary unitary class-EU The first two terms are charac- 
teristic to the chiral class. Due to this contribution, the 
number variance of the chiral unitary class iSrtwo times 
of that of the unitary class in this calculation.lLEl 
If we take the zero mode only, 



A 



{p{E 1 )p{E 2 )\ 

1 r i 

2 [4(ttE/A) 2 + (ttw/A) 2 



(p(E 1 ))(p(E 2 )) 
1 



(34) 



whichnis the asymptotics of the exact result at the ergodic 
limit .13 

For the chiral orthogonal and symplectic classes, this 
function becomes two times larger than that of the chiral 
unitary class. This is due to the cooperon contribution. 
This feature is almost the same as the traditional classes. 



III. DIAGRAMMATIC PERTURBATION 
THEORY 

In this section, we reexamine the results obtained 
above with the diagrammatic method. We are especially 
interested in the kind of diagrams that make chiral disor- 
dered systems different from traditional ones. This will 
give us some insights into the kind of quantities that are 
characteristic to the chiral classes. We work with the 
ergodic limit for simplicity. The generalization to the 
diffusive regime can be done along the same line as the 
traditional classes. Similar calculations have been done 
for systems with particle-hole symmetry] 2 ! 

We can summarize the Gaussian correlation law of the 
chiral symmetric Hamiltonian for each classes as 



(Hij Hm) 
(Hij Hm) 



Sii S k j - (S 2 )iz (S z )fcj , unitary, (35) 
Sa Skj + Sik Sij 



— (E z )jj (Ez)fcj - (Ez)ifc (£>z)lj , orthogonal, (36) 
v 2 r 

(Hij H M ) = — S k j + (cTy)ik (cTyjlj 

(E z )jy - (T, z a v ) ik (a y E z )ij , symplectic. (37) 

The effects of chiral symmetry are characterized by the 
presence of the terms including £ z . 

First, summing the rainbow-type diagrams, we can get 
the familiar form of the Green function as 



G {s) (E)) = 



tt 2 E 
8N 2 A 2 



2NA 



1- 



4^) • < 38 > 



where s = +1(— 1) for the retarded (advanced) function 
and 2N is the size of the Hamiltonian. 

We know that the quantum-interference effects are ex- 
pressed by the diffuson and cooperon diagrams. The dif- 
fuson diagram is calculated as 



n 



ij,kl ( E U E 2) 

2NA 2 



Su S. 



a okj 



1 



(^z)il (S z )fej 



'^y(G^(E 1 ))(G^')(E 2 ))\ 
Using the Green function (R8h , we have 



(39) 



/2JVAY 



V 7T / 



4iV 



■ SaS 



ilOkj 



uti7 ) (E 1 ,E 2 ) 



ij,kl 



W Suhj + i 



. (S z )ij (£z)/cj 
irE/A 

(£ z )« (S z )fej 
ttE/A 



(40) 



(41) 



n 



(+-) 

ij,kl 



(Ei,E 2 ) 



t(-,+) 



SaS, 



1 



7rw/2A 47V 



(E>z)il (^z)kj 



^7kJ , (E 1 ,E 2 ) 



ij,kl 



SaS, 



1 



7rj/2A 47V 



(42) 
(43) 

where E = (Ex + E 2 )/2 and u = E 1 - E 2 . There 
are singular contributions (1/E or 1/lu terms) and non- 
singular contributions. The remarkable property is that 
singular 1/E terms appear even in the retarded-retarded 
(advanced-advanced) channel. 

For the chiral orthogonal class, we have the cooperon 
contribution 



2iVA 2 



Sik Sij 



-(^) < {GW(E 1 ))(GM(£ 2 )) 

(^z)ik (Ez)lj 

1+( 2 ^) 2 (GM(E 1 ))(G^')(E 2 )) 



(44) 



This is the same as the diffuson contribution except the 
indices k «-> I. For the chiral symplectic class, the Pauli 
matrix a y is inserted in the cooperon contribution as 
Eq.@. 
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A. Density of states 

For the chiral unitary class, only the diffuson-type lad- 
der appears. The diagram for the lowest order correction 
of the Green function is depicted in FigJ^. The diagrams 
of this order consist of two singular-diffuson diagrams and 
the same ones with nonsingular corrections. (In Fig|l], we 
do not write the nonsingular contributions that consist 
of nine diagrams.) We get the Green function 



trG (R) (E) 



7T 

1 A 



1 



1 



(ttE/A) 2 

which coincides with the supersymmetric calculation 



(45) 



FIG. 1. Diagrams contributing to the Green function for 
the chiral unitary class. The dark-shaded regions represent a 
diffuson-type ladder. 

For the chiral orthogonal and symplectic classes, we 
can anticipate the 1/E term contribution. This is indeed 
the case and we can express the 1/E term contribution 
as Fig.^J, We note that the singular contribution comes 
from the cooperon-type ladder and the diffusion ladder 
contribution is the nonsingular one. We obtain the Green 
function as 



A 



1 



2irE/A 



16(nE/A) 2 



(46) 



We confirm that the first three terms coincide again with 
the supersymmetric calculation. 





FIG. 2. Diagrams contributing to the Green function for 
the chiral orthogonal and symplectic classes. The light-shaded 
regions represent a cooperon-type ladder. 



B. Two-point level-correlation function 

The leading order contributions to the twp-point level- 
correlation function can be expressed with the same di- 
agram as the traditional classes (two diffusons and two 



cooperons) .0 For the chiral unitary case, we take the 
diffuson-type diagram only. Calculating these diagrams, 
we can derive Eq.(CM). 



IV. CONCLUSIONS 



In this paper, we have formulated the impurity per- 
turbative expansion for chiral symmetric disordered sys- 
tems. For three chiral classes, we calculated the large- 
E and small- E behavior of the density of state and the 
large- E behavior of the two-point level-correlation func- 
tion. The deviations from the universal function are ob- 
tained in which we find the effects similar to the weak 
localization effects for the traditional classes. We calcu- 
late these deviations by considering the nonzero-A; modes 
explicitly. Alternatively, one can, in principle, integrate 
out the nonzero modes from the start and get the effec- 
tive action for the zero mode. Then, working with this 
action should give the same results. 

For the behavior at the large energy, we can interpret 
the results diagrammatically. We see that the diffuson 
and cooperon modes generate the nontrivial contribution 
to the level correlation functions. We have the diffuson 
and cooperon contributions in the advanced-advanced 
(retarded-retarded) channel in addition to the usual ones. 

For the chiral unitary and orthogonal classes we get the 
small- E behavior that could be interpreted as a weaken- 
ing of the repulsion between a pair of levels with the en- 
ergy ±E. For the chiral symplectic case, however, we get 
a result with opposite tendency whose intuitive explana- 
tion remains an open problem. These results imply the 
divergence of the density of states at the band center for 
the chiral unitary and orthogonal classes and the vanish- 
ing for the chiral symplectic class in the thermodynamic 
limit. For the chiral unitary class, this divergence-is ob- 
served by the renormalization group calculation£3 The 
calculations for the other classes are future problems. 

If once we derive the nonlinear sigma model and for- 
mulate the perturbation theory, it is easy to use this cal- 
culation to several applications such as the study of the 
wave function statistics and the renormalization group 
equations. These are currently under intensive study. 
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